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Abstract 

In previous studies we found that models with flavor-universal suppression of the neutrino-gauge 
couplings are compatible with NuTeV and Z-pole data. In this paper we expand our analysis to 
obtain constraints on flavor-dependent coupling suppression by including lepton universality data 
from W, T, IT and K decays in fits to model parameters. We find that the data are consistent with 
a variety of patterns of coupling suppression. In particular, in scenarios in which the suppression 
arises from the mixing of light neutrinos with heavy gauge singlet states (neutrissimos) , we find 
patterns of flavor-dependent coupling suppression which are also consistent with constraints from 
H 67. 
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I. INTRODUCTION 

Recent analysis of i/^ (z/^J scattering data from the NuTeV experiment at Fermilab P] 
indicates a value of the effective neutrino-quark coupling parameter g1 which deviates by 
3(7 from the Standard Model prediction (based on a global fit using non-NuTeV data). The 

n 

significance of the NuTeV result remains controversial [2], and a critical examination of the 
initial analysis is ongoing. Several groups are evaluating potential theoretical uncertainties 
arising from purely Standard Model physics which might be comparable to or larger than the 
quoted experimental uncertainty of the NuTeV result. Candidate sources of large theoreti- 
cal uncertainty include next-to-leading-order (NLO) QCD corrections j^l, NLO electroweak 
corrections j^, and parton distribution functions (especially as involves assumptions about 
sea-quark asymmetries) 0|. The effect of the first has been estimated to be comparable in 
size to the NuTeV experimental uncertainty, while the latter two might give rise to effects 
comparable in size to the full NuTeV discrepancy with the Standard Model. Elucidation 
of the actual impact of these effects on the NuTeV result awaits a reanalysis of the NuTeV 
data. However, it remains a distinct possibility that the discrepancy with the Standard 
Model prediction is genuine and that its resolution lies in physics beyond the Standard 
Model p. It is this possibility that we investigate here. 

In a previous paper [3], we demonstrated that the Z-pole data from e+e~ colliders ^] and 
the z/^ (z/^) scattering data from NuTeV [l| are compatible if (1) the Higgs boson is heavy 
and (2) the ZuiUe and Wiue {i = e, fi, r) couplings are suppressed by a factors (1 — ee) and 
(1 — ee/2), respectively. We also showed that such suppressions could arise from neutrinos 
mixing with heavy gauge singlet (neutrissimo) states 9l. llOl IllL Il2|. 

In Ref. ^, it was assumed that the suppression parameters were flavor universal: Se = 
e^ = Sr = e. The value of e required to fit the data was 

£ = 0.0030 ±0.0010 . (1) 

However, in seesaw models of neutrino masses and mixings such a large universal e 
implies a prohibitively large rate of /i — >^ 67 Q| . To bring the models into agreement 

experiment the assumption of universality must be relaxed: either or but not both, 
must be strongly suppressed Further, in most models flavor-universal suppressions 

require considerable fine tuning. It is thus natural to ask what patterns of flavor non- 
universal suppressions are consistent with the data. If the suppression parameters can be 



flavor-dependent, one must also ask whether the preferred values of the are all positive, 
i.e. are all the neutrino-gauge couplings suppressed? Negative indicates an enhancement 
of the WCvi and Zvfyi couplings which cannot be arranged via neutrino mixing. 

In addition to the Z-pole and NuTeV data, there is a wealth of experimental data bound- 
ing lepton universality violation in the charged channel from VT, vr, i^, and r-decays jl^j- In 
the following, we analyze the constraints that Z-pole and NuTeV data and the lepton univer- 
sality bounds impose on neutrino-mixing models by fitting the data with flavor- dependent 
suppression parameters Si [i = e, /i, r) along with the S, T, and U oblique correction pa- 
rameters [1^. We perform fits in which all six parameters float independently, and we also 
fit to models in which one or more of the Ei are assumed to be strongly suppressed. As in 
the flavor-universal case, the data require a negative T parameter and a positive U parame- 
ter. However, we find that the data are consistent with a variety of patterns of suppression 
parameters, including patterns compatible with /i — 67 data. 

II. CONSTRAINTS ON LEPTON UNIVERSALITY 

Here, we survey current experimental constraints on lepton universality. For a compre- 
hensive review on the subject, see Ref. We parametrize the couplings of the VT^'s with 
the leptons as 



The Standard Model assumes Qe = Ofi = Or = 9- Although experimental limits on the 
ratios g^/Qe, gr/gii, and gr/ge have been calculated and tabulated as recently as fall 2002 
by Pich we repeat the exercise here to incorporate more recent data and to obtain the 
correlations among the limits necessary for our analysis. 

A. VF- decay 
The decay width of the W at tree level is 
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The branching fractions of the W into the three lepton generations have been measured at 
LEP II to be (Ref. 3, page 74) 



B{W ^ 

B{W TVr) 



10.59 ±0.17% 
10.55 ±0.16% 
11.20 ±0.22% 



(4) 



with correlations shown in Table HI From this data, we find (Ref. p| , page 73) 



B{W^ fiiyf,)/B{W ^ ez7e) 
B{W TVr)/B{W eVe) 



0.997 ±0.021 , 
1.058 ±0.029 , 



(5) 



with a correlation of ±0.44 between the two ratios. Using Eq. Q, this translates into 

{9j9e)w = 0.999 ± 0.011 , 

igr/ge)w = 1.029 ± 0.014 , (6) 

with a correlation of ±0.44. The central values have shifted slightly from Ref. 0| due to 
the update of the W branching fractions from LEP II. 





B{W ^ eve) B{W ^ 


B{W TUr) 


B{W - 




1.000 0.092 


-0.196 


B{W - 




1.000 


-0.148 


B{W - 






1.000 



TABLE I: Correlations among the W branching fractions measured at LEP II (Ref. [a], page 182) 



B. r and fi decay 

The decay widths of the r and jj, into lighter leptons, including radiative corrections |3, ll 
are: 
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phase space 


W propagator 


photon 


r(T li.i'^h'T) 


f {ml/ml) = 0.9726 


6^ = 1.0003 


6:^ = 0.9957 


r(T ^ eVe Vr) 


f {ml/ml) = 1.0000 








f (ml/ml) = 0.9998 


S'^ = 1.0000 


5/; = 0.9958 



TABLE II: The corrections to the leptonic decay widths of the r and /x. 
in which f{x) is the phase space factor 

f{x) = l-8x + 8x^ -x"^ - 12x^ In x , 
6^ is the W propagator correction 

6^ is the radiative correction from photons 

^ ^ ajme) /25 

and the values of the running QED couphng constant at relevant energies are 



(8) 



(9) 



TX 



(10) 



mj = In —^i + — ^ 136.0 

Stt OTT 



a ^imT-) ~ 133.3 



(11) 



The numerical values of these corrections are shown in Table ITTl The ratios of the coupling 
constants can be extracted using the relations 



r(r 

Tir- 



■ eVeVr (7)) 
/^^^M^^r (7)) 



B{t ^ eVeVr (7)) 



B{t ^ fXU^Ur (7)) 



ge f {ml/ml) ' 
^ gl ml f (ml/ml) 5; 

9e /("^e/K) ^1 



(12) 



r(/i ^ ez/eZ/^(7)) 

The latest world averages for the quantities appearing in these equations are listed in Ta- 
ble Hill which yield 



{9T/9e)Tfi 



0.9999 ±0.0021 , 
1.0004 ±0.0022 , 



(13) 



with a correlation of 0.51 due to the inputs B{t — > fiu^Uj. (7)) and common to both 
ratios. 
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Observable 


World Average 


Reference 


me (MeV) 


0.510998918 ± 0.000000044 


[21] 


(MeV) 


105.6583692 ± 0.0000094 


[2JJ 


(s) 


(2.19703 ±0.00004) x 10"^ 


[21] 


nir (MeV) 


1776.99l|]:26 


[2JJ 


Tr (s) 


(290.6 ±0.9) X 10-15 


[23] Figure 1 


B{t — > e z/g i^t) 


17.823 ±0.051% 


[24] Figure 8 


fi(r llU^Ur) 


17.331 ±0.054% 


[24] Figure 9 


B{t — > vr l^r) 


10.975 ± 0.065 % 


[21], [23] Figure 3, [24] Table 3 


B{t KVr) 


0.686 ±0.023% 




rrin (MeV) 


139.57018 ± 0.00035 


[21J 


Tn (s) 


(2.6033 ± 0.0005) x 10"^ 




5(7r Hh'^) 


99.98770 ±0.00004% 


[2JJ 


B{7r eVe) 


(1.230 ±0.004) X 10"^ 


[21] 


niK (MeV) 


493.677 ± 0.016 


[2JJ 


TK (s) 


(1.2384 ±0.0024) x 10"*^ 


[21J 


B{K^^l9^) 


63.43 ±0.17% 


[21] 



TABLE III: The world averages of masses, life times, and branching fractions used in this analysis. 
The branching fractions subsume the decays with 7's. 

C. Pion and r decay 



At tree level, the widths of charged vr-decay into leptons are 



r(7r ez/e) 

r(vr ^ ^iz>^) 
while that of r-decay into tt Vt is 



2567r " " V 

T- ^.Mtt 1 T 

2567r " m2 y 



^rfl'nd /tt 3 ( i 
5127r ^ 1^ 



m: 



(14) 



(15) 



where Qud = and the pion decay constant f-,, is normalized as (Ref. [20], page 439) 

(OlM7^75ci(0)|7r'(q))=^g^A. (16) 
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Processes 


Constraint 


W - 




{gti/de 


w = 


0.999 ±0.011 


w - 


— > U Vn 


\HT 1 He / 


VV 


1.029 ± 0.014 


w - 


T Ur 








A* - 




{9tJ./9e] 


r = 


0.9999 ±0.0021 


T — 


> eUeVj- 


\HT ! i/e / 




1.0004 ± 0.0022 


r — 










TT — 




(Qii /Op' 


TT = 


1.0021 ± 0.0016 


IT — 






'ttt — 


1.0030 ± 0.0034 


T — 


> TT Vr 








K - 






)Kt = 


= 0.979 ± 0.017 


T — 











TABLE IV: Limits on lepton universality from various processes. 



Taking ratios, we find 

T{n ^ eue) _ B{TT^eVe) _ gl ml {I - ml/ml)'^ 

9l ml [l-ml/mlY ' 



V{t^ttv^) _ B{t ^ TTU^) _ gl ml {l-ml/m^f 



TT 



gl 2mlm^ (1 -ml/miy 



(17) 



Radiative corrections to these relations have been calculated in Ref. j^l and modify them 
to 



Rt/-k 



with 



5R, 



B{n ^ eue (7)) _ 
fi(7r ^/iz/^(7)) 
B{t ^nur (7)) 

Tr 5(7r ^/iZ/^(7)) 



-0.0374 ±0.0001 



Re/,1 (1 + ^Re/fi 
= Rt/it (1 + ^Rt/tt 



SRrU = ±0.0016+°-°°°9 



(18) 



(19) 



t/tt — -ru.uuxu_o,ooi4 

The uncertainty in these corrections is due to the uncertainty from strong interaction effects. 
With these relations and the experimental data listed in Table IIIH we obtain 



{gfi/ge)7r 
(fi'r/fi'M)7r7 



1.0021 ±0.0016 , 
1.0030 ± 0.0034 . 



(20) 
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The correlation between the two is virtually zero due to the accuracy of the common inputs 
rrifj,, mTT, and B^n /i^'^)- There is a correlation of +0.33 between ((7^/(7^)^^ and {gT/ge)Tfi 
of Eq. (|13|) arising from the common inputs and m^. A few comments are in order: 

• Our limit on {gfi/ge)n differs from that of Pich and Silva 3| who use for the value of 
B{n — > eUeil)) the weighted average of the resuhs from TRIUMF H and PSI Q, 
(1.2310 ± 0.0037) X 10~^. If we use this value instead of the average from the Review 
of Particle Properties J2l| listed in Table mH we obtain {g^/ ge)-K = 1.0017 ± 0.0015 in 
agreement with Ref. j26[. 

• The experimental value of i?(r thO listed in Table ITTTl is the average of CLEO and 
the four LEP experiments. CLEO Q, OPAL 0, DELPHI and L3 Q report 
the semi-exclusive branching fraction B{r —>■ hv^-), where h = ti 01 K, a,s 

CLEO : B{t hv^) = 11.52 ± 0.05 ± 0.12 % 

OPAL : = 11.98 ± 0.13 ± 0.16 % , , 

(21) 

DELPHI : = 11.601 ± 0.120 ± 0.116 % 

L3 : = 12.09 ± 0.12 ± 0.10 % 

The CLEO and OPAL values are published and used in the average of the Review of 
Particle Properties j2l|. Adding the statistical and systematic errors in quadrature 
and taking the weighted average of these four numbers, we obtain 

average without ALEPH: B{t ^ hi^r) = 11-752 ± 0.079 % . (22) 



As noted by Can [2J|, the agreement among these four measurements is poor 
(x^/d.o.f. = 9.9/3, where d.o.f. is degrees of freedom). This stands in stark con- 
strast to the situation in fall 2002 when the agreement among CLEO and the four 
LEP experiments was much better (x^/d.o.f. = 5.09/4) 3^. The source of the differ- 
ence is the new L3 value 32 1 which has a much higher central value and smaller error 

n 

3ar than before [3J|. Subtracting the world average B{t Kv^-) = 0.686 ± 0.023% 



2l|, we obtain 

average without ALEPH: B{t nUr) = 11.066 ± 0.082 % . (23) 
ALEPH (js^. Table 3) reports the value of the exclusive branching fraction B(t 



ttUt-i as 



ALEPH 



B{t TiUr) = 10.828 ± 0.070 ± 0.078 



(24) 





{g^ilge)w {gT/ge)w 


{g^/ge)T {9Tlge)T,ji {gf,/ge)-K 


{gTlg^)-KT 


{gT|g^l)KT 


{9^l|ge)w 


1.00 0.44 








{gr/ge)w 


1.00 








{gfil9e)T 




1.00 0.51 0.00 


0.00 


0.00 


{gT/ge)Tti 




1.00 0.00 


0.33 


0.07 


{gfilge)-K 




1.00 


0.00 


0.00 


{gT/g^i)TTT 






1.00 


0.04 


{gT/g^i)KT 








1.00 



TABLE V: Correlations among the lepton universality constraints. 



which does not agree particularly well with Eq. ()23p either. Although this ALEPH 
value is excluded from the world average by Gan in Ref. Q] as preliminan^, we include 
it in our analysis since it was included in the previous analysis by Pich jl2| (with the 
caveat that it is subject to change). The weighted average with Eq. is 

world average: B{t txu^) = 10.975 ± 0.065 % , (25) 

which is the value used to obtain Eq. (pUj) . The associated x^/d-o.f. is 13.1/4, so is 
unimproved with the inclusion of the ALEPH result. If we exclude the ALEPH value 
and use Eq. ^ instead, we obtain {g^/g^)^^ = 1.0072 ± 0.0041. 

The current state of agreement among the data determining B{t ^ ttz/,-) is clearly 
unsatisfactory. Additional data, perhaps from new experiments at CLEO j^^, are 
needed to provide a definitive value. 



D. Kaon and r decay 

Paralleling the treatment of pion decays, we can extract gr/ Q^i from kaon decays. The tree 
level decay widths involving kaons are 

.2 \ 2 

" 



r(/f-A.^,.) = li?#™>Wi-3^) . (26) 

and 



2567r - V m-j^ 



2„2 

" (27) 



r(r ^ z/.) = ^ -^mM 1 - ^4 



5127r M ml 
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where Qus = glVusl, and the kaon decay constant Jk is normahzed as (Ref. j^], page 439) 

(0|M7M75s(0)|ir-(q))=zg^/^. (28) 

Taking the ratio yields 



r(r ^ Kur) _ TK BjT^Kur) _ gl ml (1 - mf^/m^)^ 
r(7r /i z/^) B{K i/^) 2mlmK (1 - ml/m\f 



Radiative corrections modify this to 



with 



-Kr/X - ^TTT / - K/K + Onr/K) , (.3UJ 



6Rr/K = +0.0090^°:°°^^ . (31) 
Using this relation and the data listed in Table ITTTl we obtain 

(W^7/.)i^r = 0.979 ±0.017, (32) 

which agrees with Ref. [l^]. This has a correlation of +0.07 with [g^/ g,,)^^ of Eq. (fT^ . and 
a correlation of +0.04 with (g^/g^)^^ of Eq. arising from the common inputs Tr and 

We tabulate our results in Tables HYI and El 



III. Z-POLE, NUTEV, AND W MASS DATA 

For the Z-pole and NuTeV data, we use the same set as in Ref. 0, namely Fiept, Tmv/^\ept, 
and sin^ 6'^'^* from e~^e~ colliders, g1 and g]^ from NuTeV. Of these, only the value of sin^ 6^^^^ 
has been updated since our last analysis in Ref. The W mass has also been updated 
by LEP-II. We list the values in Table IVTl There is a correlation of 0.17 between Fiept and 
Einv/Eiept; other correlations are negligible. 



IV. THE CORRECTIONS 



Suppression of the neutrino-gauge couplings modifies the relation between the Fermi 
constant Gp and the muon decay constant to 

Gf = G,{i + . (33) 
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Observable 



Mz 

Tlept 

-Tinv/ricpt 
-.2 /jlcpt 



sm 



9l 
Mw 



SM prediction 


Measured Value 


Reference 


Input 


91.1875 ± 0.0021 GeV 


fa,] page 8 


84.034 MeV 


83.984 ± 0.086 MeV 


[8] page 9 


5.972 


5.942 ±0.016 


[8] page 8 


0.23133 


0.23150 ±0.00016 


[8] page 142 


0.3040 


0.3002 ±0.0012 


Average of \\] and \22] 


0.0304 


0.0310 ±0.0010 


Average of W] and [22j 


80.399 GeV 


80.426 ± 0.034 GeV 


[8] page 146 



TABLE VI: The observables used in this analysis in addition to the lepton universality data. The 
measured value of sin^ 9^^^^ and the W mass have been updated in Ref. 8] since the analysis in 
Ref. 7]. The SM predictions are ZFITTER 36] outputs with inputs of A/top = 178.0 GeV 
Mffiggs = 115 GeV, as{Mz) = 0.119, and Aaj^^^d = 0.02755 Q- 



Since is used as an input in calculating the SM predictions, all observables whose SM 
prediction depends on Gp will receive this correction through Gp- Of the observables that 
measure the ZviVg and Wiv(, vertices directly, the Z invisible width is corrected by an 
additional factor of 

2 

1 - 3 (^e + + ^r) , (34) 

while the NuTeV parameters g\ and g\ receive an additional correction of (1 — e^). The 
dependence of the observables on the oblique correction parameters 5, T, and V can be 
found elsewhere jisl ]. 

Numerically, the observables are corrected as follows: 

Tlept 



[riept]sM 
Tinv/riept 



[rinv/riept]sM 

[sin^ ^if *]sM 



9l 



lJsm 
9\ 



[fi'^lsM 



w 



[M, 



1 - 0.0021 S ± 0.0093 T ± 0.60 ± 0.60 e^, , 
1 ± 0.0021 S - 0.0015 T - 0.76 £e - 0.76 - 0.67^^ 
1 ± 0.016 S - 0.011 T - 0.72 Se - 0.72 e, 
1 - 0.0090 5 + 0.022 T + 0.41 - 0.59 
1 + 0.031 S - 0.0067r - 1.4£e - 2.4^^ 
1 - 0.0036 S ± 0.0056 T ± 0.0042 f/ ± 0.11 + 0.11 



'At ; 



SM 



(35) 
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Here, [*]sm is the usual SM prediction of the observable * using as input. 

Despite the fact that six observables are available to the fit, this set of data is not sufficient 
to fix all six parameters. This is because the ratio g^/gpi depends on the fit parameters only 
through sin^6''^*, and thus can only constrain the exact same linear combination of S, T, 
Ee, and as sin^^^'^*. In fitting the parameters to the observables, the linear combination 

aT + 2{e^~2ee) , (36) 

remains unconstrained. Therefore, we can constrain only five of the six fit parameters with 
the Z-pole and NuTeV data. 

The linear combinations constrained by the lepton universality bounds from W, r, vr, 
and i^'-decays are 



5V 


= 1 + 


£e 




9e 




2 




= 1 + 


£m 


- Sr 








2 




= 1 + 




- St 


9e 






2 



(37) 



Since there are only two independent observables but three fit parameters, only two inde- 
pendent linear combinations can be simultaneously constrained by the lepton universality 
data. However, when these bounds are combined with the Z-pole and NuTeV data, all three 
El can be constrained independently. 

V. FITS TO LEPTON UNIVERSALITY CONSTRAINTS 

The three ratios of the coupling constants contain only two independent degrees of free- 
dom, since the third is the product of the first two. Equivalently, when parameterizing the 
ratios in terms of the differences of the e^, as in Eq. ()37|1 . only two [any two) sets of the 
differences can be taken as free parameters; the third can be expressed as the difference 
between them. We can combine the seven pieces of experimental data in Table HVl bv fitting 
with any two of the three parameters 





= Be 


-Bf, = 


AgT -|- ^fJ,T 5 




A^, 


= 


-Er = 


AgT Ag^ , 




Ae. 


= Se 


-Er = 


Ae/j ^fiT ■ 


(38) 
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Ae^ 


A^, 


Aer 




1.00 


-0.27 


0.37 


A^^ 




1.00 


0.80 


Aer 






1.00 



TABLE VII: Correlations among the A's from fit. 

We obtain 

Ae^ = 0.0022 ±0.0025 , 
A^^ = 0.0017 ±0.0038 , 

Ae^ = 0.0039 ± 0.0040 , (39) 

with correlations shown in Table IVlTl In terms of the coupling constant ratios, this translates 
to 

{gjg,) = 1.0011 ±0.0012 

{g^/g^) = 1.0019 ±0.0020 (40) 

with a correlation of 0.37. The quality of the fit is unimpressive: the is 8.4 for (7- 
2)=5 degrees of freedom. The largest contribution is from {gr/ ge)w which contributes 4.6. 
The region of AeT--A^^ parameter space preferred by the fit is shown in Fig. ^ The 90% 
confidence contour preferred by the VT-decay data hardly overlaps with that of the r-decay 
data, which causes the large x^. 

Since the objective of this paper is to determine whether the Z-pole and NuTeV data 
are compatible with lepton non-universality, it is problematic that the lepton universality 
constraints are not clearly compatible among themselves. The set of coupling ratios we 
consider here has an intrinsic of 8.4 (or 10.8, if the ALEPH value is excluded from the 
calculation of -B(r ^ t^Vt)) which cannot be mitigated in our model. This is in addition to 
the large associated with B{t t^Vt)-, discussed previously. Further experiments may 
provide the ultimate resolution of the tension in the data. For now, to prevent this large 
among the lepton universality data from obscuring their compatibility with the Z-pole and 
NuTeV data, we will use the average values obtained in Eq. pn|) in our subsequent fits with 
the caveat that the pair hides a large x^. 
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VI. FITS TO Z-POLE, NUTEV, AND LEPTON UNIVERSALITY DATA 



We fit the expressions in section IV to the Z-pole, NuTeV, and W mass data listed in 
Table EU and the lepton universality constraint Eq. (j4(J|) . The S, T, U parameters were 
used in all fits. Of the three we performed fits with the following eight combinations of 
fit parameters: 

A. fit with a flavor- universal e {e^ = = Er = e), 

B. fit with all three parameters Se, s^, and e,-, 

C. fit with Be and e^, 

D. fit with £e and Sr, 

E. fit with and Sr-, 

F. fit with Be only, 

G. fit with £^ only, 

H. fit with £r only. 

Fit A with flavor-universal e is the one performed in Ref. [3| (without the lepton universality 
constraints). We include it here as a benchmark against which to compare the fiayor- 
dependent fits. The reference Standard Model values were calculated using ZFITTER 
with the inputs Mz = 91.1875 GeV 
0.119, and Aa'^^iMz) = 0.02755 



Mh = 115 GeV, Mt = 178.0 GeV |33], a,(Mz) 



The results of these fits have been tabulated in Table IVIIH with correlations among the 
fit parameters shown in Table |X| As the values of in Table IVIIII indicate, the quality of 
the fits A, B, C, D, and F is excellent, while fits E and G are only marginal and fit H fails. 
The largest contribution to the overall fo^^ fits E and G is from {g^/ge) (5.1 for E and 
5.3 for G) which indicates that these fits are not compatible with lepton universality. For 
fit H, the largest contributions to the overall is from the NuTeV observable g1 (5.5) and 
[gr/ge) (3.0) which indicates that neither NuTeV nor lepton universality are accommodated. 
Comparisons of fits B and C, D and F, E and G show that including er in the fits does little 
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fit parameters 


A 


B 


C 


D 


S 


-0.01 ±0.10 


0.00 ±0.10 


0.00 ±0.10 


-0.04 ±0.10 


T 


-0.48 ±0.15 


-0.56 ±0.16 


-0.56 ±0.16 


-0.45 ±0.15 


U 


0.55 ±0.16 


0.62 ±0.17 


0.62 ±0.17 


0.51 ±0.16 




0.0030 ± 0.0010 


0.0048 ± 0.0018 


0.0049 ± 0.0018 


0.0050 ± 0.0018 




0.0030 ± 0.0010 


0.0027 ± 0.0014 


0.0027 ±0.0014 


— 




0.0030 ± 0.0010 


0.0007 ± 0.0028 





0.0012 ±0.0028 




2.4 


0.91 


0.97 


4.4 


d.o.f. 


4 


2 


3 


3 


xVd.o.f. 


0.6 


0.5 


0.3 


1.5 


large 











fit parameters 


E 


F 


G 


H 


S 


-0.03 ±0.10 


-0.04 ±0.10 


-0.03 ±0.10 


-0.08 ±0.10 


T 


-0.30 ±0.13 


-0.46 ±0.15 


-0.30 ±0.13 


-0.18 ±0.12 


U 


0.38 ±0.15 


0.52 ±0.16 


0.37 ±0.15 


0.25 ±0.13 






0.0051 ± 0.0018 








0.0028 ± 0.0014 




0.0029 ± 0.0014 






0.0021 ± 0.0027 






0.0026 ± 0.0027 




7.8 


4.5 


8.3 


11.6 


d.o.f. 


3 


4 


4 


4 


xVd.o.f. 


2.6 


1.1 


2.1 


2.9 


large 


5.1 from {g^/ge) 
2.5 from {gr/ge) 




5.3 from (fif/i/Se) 


5.5 from g\ 
3.0 from {gr/ge) 



TABLE VIII: The results of the fits. 



to improve the overall x^. Indeed, the per degree of freedom in actually worse with the 
inclusion of e^. 

The constraints placed on the fit parameters by each observable are illustrated in figure El 
The la bands in each 2 dimensional plane are plotted assuming that all other fit parameters 
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fit parameters 


E' 


F' 


S 


-0.03 ±0.10 


0.00 ±0.10 


T 


-0.50 ±0.17 


-0.38 ±0.14 


U 


0.56 ±0.18 


0.45 ±0.15 


£e 


0.0058 ± 0.0023 







— 


0.0048 ± 0.0017 




4.7 


2.9 


d.o.f. 


4 


4 


xVd.o.f. 


1.2 


0.74 



TABLE IX: The results of tlie fits using tlie lepton universality constraint from ^ and r-decay only. 

are set to zero. The gray elUpses are the 68% and 90% confidence regions for fit C, i.e. 
the five parameter fit with T, U, Ee, and e^, projected to each plane from the full five- 
dimensional parameter space. In the case of the Mw band, since Mw serves only to fix U, 
it exerts no statistical 'pull' on the other fit parameters; also, since U is fixed by a single 
observable we have not included figures with U on an axis. We have also omitted projections 
onto planes involving since serves little in improving the quality of the fits and since, 
other than the lepton universality constraint on (fi'r/fl'e), it is constrained by only Ti^^/Tiept- 
Figure El clarifies the reason for the failure of fit H. From figureEfa), we see that the NuTeV 
observable gj^ prefers a negative T. To maintain the agreement between the SM predictions 
and the Z-pole data, the effect of negative T must be absorbed by a corresponding shift 
in Gp, Eq. (jH^ . by making Ee and/or positive (as indicated in figures Efd) andEfe)). 
However, since Eg and are both constrained to zero in fit H, it cannot accommodate gj^. 
Further, the measured value of Finv/Ficpt is smaller than the SM prediction, which demands 
positive Er, while (gr/ge), Eq. (jlOI), prefers negative Er. Thus, H cannot satisfy {gr/ge) 
either. 

For fits E and G, in which Ee is constrained to zero, the effect of a negative T is absorbed 
by a positive e^. However, the experimental value of (g^/ge) prefers e^ negative. A tension 
thus remains between the electroweak data and the lepton universality data in these fits. 
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VII. DISCUSSION AND CONCLUSIONS 



The electroweak data are well-fit by several of the patterns of neutrino-gauge coupling 
suppressions considered. In all cases considered, the best-fit values of the are positive, i.e. 
neutrino-gauge couplings are suppressed with respect to the Standard Model, as is required 
in models of neutrino mixing. In models in which Se is allowed to be non-zero (A-D,F) the 
fit to the data is good, and the fit improves if is allowed to be non-zero as well (B,C). 
The fit quality is degraded for models in which Se is set to zero (E,G,H), and the fit with Sr 
alone (H) is poor. In general the overall is insensitive to the presence of as a degree of 
freedom in the fit. The data prefer the model with only Sf, non-zero (F) to the model with 
only non-zero (G). 

Since the _/i > 67 data from MEGA 3^ demands that either or is strongly sup- 



pressed [10|, [111, 1121, 7^ 0, ~ seems to be the solution preferred by current data. 



However, we stress that the inconsistency within the lepton universality data makes any 
such conclusion tentative. For example, fits using only the lepton universality constraint 
from yU and r-decay, Eq. (jl3j) . which is free of QCD uncertainties, indicate that the fit with 
only is superior to a fit with only e^, as shown in Table \iW Therefore, future improve- 
ments in the lepton universality data {e.g. better determination of the r lifetime by Belle 



and Babar 



23|, measurement of B{tx —>■ eVe) at the 0.2% level by PIBETA [40], etc.) may 



ultimately provide a different conclusion. In figure Eff), the current 90% contour overlaps 
with the Be axis but not with the axis. If the region preferred by the lepton universality 
data (dashed contour) is shifted toward the axis, eg ~ 0, 7^ may become a viable 
solution also. 

Langacker 0] has noted that the observed violation of unitarity in the CKM matrix 42 1 



will be aggravated by suppressions of neutrino-gauge couplings. However, if the suppression 
parameters are allowed to break universality, it is only a non-zero that aggravates the 
CKM unitarity problem. Thus the CKM unitary data actually prefers the 7^ 0, ~ 
solution (in the sense that it does not make the problem worse). An improved determination 
of \Vud\ is expected from the UCN-A experiment at LANL in the near future 4^. 

The fits A, B, C, D, and F with excellent x^'s require T to be negative by more than Scr, 
U to be positive by more than 3cr, while 5* is within la of zero. As discussed in Ref. j?!, the 
5* and T parameters can be accommodated within the Standard Model by increasing the 
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Higgs mass to several hundred GeV. The large U parameter arises in part from discrepancy 
between the Standard Model prediction for the W mass and in part from the shift due to the 
other fit parameters. Neutrino mixing alone does not account for this discrepancy between 
the predicted and observed values of the W mass; the U parameter appears to require new 
physics. Whether a large U parameter can be generated without a correspondingly large T 
parameter in some model is an open question that needs to be addressed. 



44 



47 



The constraint on the suppression parameters {i = e, /i, r) from muon g — 2 
is weak 46|. Further constraints may be obtained from u to e conversion in nuclei 
and _ant™n,un, „.cmation T..e wU, be d.cu.ed . a fu.u. 

publication. 
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TABLE X: Correlations among the fit parameters for fits A thought H. 
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FIG. 1: The limits on Agr and A^.,- from (a) T^-decay, (b) r-decay, (c) tt and K-decay, and (d) all 
decays combined. The la bands are shown for each coupling constant ratio ignoring correlations. 
The shaded areas represent the 68% (dark gray) and 90% (light gray) confidence contours including 
correlations. 
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S ee 

FIG. 2: The 68% and 90% confidence contours projected onto various planes for the five-paremeter 
fit with S, T, U, Ee and e^. The bands associated with each observable show the la limits in the 
respective planes. The origin is the reference SM with Mtop = 178.0 GeV and Mniggs = 115 GeV. 
The dashed contour in 2(f) is the 90% lepton universality bound. 
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